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Abstract 



We study the distribution and various properties of exponential functionals of hypergeometric 
Levy processes. We derive an explicit formula for the Mellin transform of the exponential functional 
and give both convergent and asymptotic series expansions of its probability density function. As 
applications we present a new proof of some of the results on the density of the supremum of a stable 
process, which were recently obtained in [25] and [23]. We also derive some new results related to 
(i) the entrance law of the stable process conditioned to stay positive, (ii) the entrance law of the 
excursion measure of the stable process reflected at its past infimum and (iii) the entrance law and 
the last passage time of the radial part of n-dimensional symmetric stable process. 
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1 Introduction 

Exponential functionals of Levy processes play a very important role in various domains of probability 
theory, such as self-similar Markov processes, random processes in random environment, fragmentation 
processes, branching processes, mathematical finance, to name but a few. In general, the distribution of 
the exponential functional of a Levy process X = {Xt,t > 0) with lifetime (, defined as 



can be rather complicated. Nonetheless, this distribution is known explicitly for the case when X 
is either: a standard Poisson processes, Brownian motion with drift, a particular class of spectrally 
negative Lamperti stable (see for instance [18, 28, 34]), spectrally positive Levy processes satisfying the 
Cramer condition (see for instance [36]). In the class of Levy processes with double-sided jumps the 
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distibution of the exponential functional is known in closed form only in the case of Levy processes with 
hyperexponential jumps, see the recent paper by Cai and Kou [14]. An overview of this topic can be 
found in Bertoin and Yor [7]. 

For many applications, it is enough to have estimates of P(/ < t) as t — ?■ 0^ and P(/ > t) as t — ?■ +oo. 
However it is quite difficult to obtain such estimates in the general case. The behaviour of the tail 
P(J > t) has been studied in a general setting (see for instance [18, 19, 22, 32, 37, 38]). On the other 
hand, the behaviour of P(/ < t) has been studied only in two particular cases: when X has exponential 
moments and its Laplace exponent is regularly varying at infinity with index 6 G (1,2) (see [33]) and 
when X is a subordinator whose Laplace exponent is regularly varying at (see [9]). 

At the same time, the problem of finding the distribution of the supremum of a stable process has 
also stimulated a lot of research. The explcit expressions for the Wiener-Hopf factors for a dense set 
of parameters were first obtained by Doney [21]. In the spectrally positive case a convergent series 
representation for the density of supremum was first obtained by Bernyk, Dalang and Peskir [4], and 
a complete asymptotic expantion was derived by Patie [35]. The general case was treated recently in 
[25] and [23], where the authors have derived explicit formulas for the Wiener-Hopf factors, the Mellin 
transform of supremum and also convergent and asymptotic series representations for the density of the 
supremum. 

In this paper we pursue two goals. First, we will study exponential functionals of hypergeometric 
processes, we will obtain the Mellin transform and both convergent and asymptotic series representations 
for the density of the exponential functional. This gives us the first explicit results on the exponential 
functionals of Levy processes which have double-sided jumps of infinite activity or infinite variation. 
The technique that we develop to obtain these results may be of independent interest, as it provides a 
very simple and straightforward method to derive well-known results on the exponential functionals of 
Brownian motion with drift and the less-known (but very interesting) recent results by Cai and Kou [14] 
for processes with double-sided hyper-exponential jumps. 

The second goal is to prove that hypergeometric processes include the Lamperti-stable processes, 
which will allow us to prove interesting results on fiuctuations of stable processes with the help of the 
Lamperti transformation. In particular, we will obtain a new proof of some of the results in [25] and [23], 
which is more straightforward and less technical. We will also derive some new formulas related to the 
density of the entrance law of the excursion measure of the stable process refiected at its past infimum 
and will obtain some several results related to the n-dimensional symmetric stable process. 

The paper is organized as follows: in Section 2 we introduce hypergeometric processes and establish 
the connections between this class and the Lamperti-stable processes. In Section 3 we study the Mellin 
transform of the exponential functional of the hypergeometric processes and in Section 4 we derive the 
convergent and asymptotic series representations for the density of the exponential functional. Finally, 
in Section 5 we present some applications of these results to fiuctuations of stable processes. 

2 Hypergeometric and Lamperti-stable processes 

Hypergeometric processes were first introduced in [29] and, more generally, in [26]. These processes were 
originally constructed using Vigon's theory of philantropy (see [40]) and they provide examples of Levy 
processes with an explicit Wiener-Hopf factorization. The class of processes which we will study in this 
paper should be considered as a generalization of the hypergeometric processes studied in [29] and of the 
Lamperti-stable processes, which were introduced by Caballero and Chaumont [10]. 
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We start by defining a function ip{z) as 

r{i-p + ',-z) T{0 + 7 + 2) 
"^'^' = - r(i-^-.) r(/j + .) ■ 

wliere 7, /3, 7) belong to the admissible set of parameters 

^ = {/3 < 1, 7 e (0, 1), /3 > 0, 7 e (0, 1)}. (2) 

Our first goal is to show that tp{z) is the Laplace exponent of a (possibly killed) Levy process X, that is 
ip{z) = lnE[exp(2;Xi)]. This process and its properties will be presented in the next Proposition. From 
now on we will use the following notation 

r/ = 1 - /3 + 7 + /3 + 7. (3) 

Proposition 1. 

(i) The function ip{z) defined by (1) is the Laplace exponent of a Levy process X. The density of its 
Levy measure is given by 



^^"^^ e-(^-^+^)"2Fi (1 + 7, ^; ^ - 7; e-") , if X > 0, 



TT X = < 



r(r/-7)r(-7) 

r(r^) 



(4) 



g(/3+7)x^^^ (1 + ^, r/; r/ - 7; e") , if x < 0, 



I r(r/-7)r(-7) 
where 2F1 is the Gauss hypergeometric function, 
(a) When (3 <1 and /3 > the process X is killed at the rate 

' r(i-fl) r(/3) 

When (3 = 1 and /3 > {/3 < 1 and /3 = 0} X drifts to +00 {— C)o} and 

Eix.i = imm. /e|.y,i = _r(7)r(i-^ + 7) 
' ' r(f)) \ ' " r(i-f)) 

When 13 = 1 and /3 = the process X oscillates. 

(Hi) The process X has no Gaussian component. When 7 + 7 < 1 {1 < 7 + 7 < 2} the process has paths 
of bounded variation and no linear drift {paths of unbounded variation} . 

(iv) We have the Wiener-Hopf factorization —ip{z) = K,{q, —z)k{q, z) where 
(v) The process X = —X is a hypergeometric process with parameters (1 — /3, 7, 1 — (3, 7) . 
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Proof. First let us prove (i). Let X^^^ be a general hypergeometric Levy process (see section 3.2 in [26]) 
with parameters 

1 1 

a = d = ki = 6i = 62 = 0, /3 = 1, Ci = — — r, C2 



tti = /3, ^2 = 1 - /3, 7i = 7, 72 = 7, h 



r(-7)' r(-72) 
r(/3 + 7) 
r(/3) 



This process is constructed using Vigon's theory of philantropy (see [40]) from two subordinators H^^'^ 
and H^^^ which have Laplace exponents 

K(^)(g, z) = K{q, z) - K(g, 0), k^^\q, z) = k{q, z) 

where ii{q,z) and k{q,z) are given by (5). We see that the Lapalace exponent ip^^\z) of process X^^^ 
satisfies 

^(i)(^) = ip(z) + kk{q,z) 

where we have denoted k = K{q,0). Therefore the process ^ has the same distribution as Xt — H^t, 
in particular the distribution of positive jumps of X^^^ coincides with the distribution of positive jumps 
of X. From [26] we find that the Levy measure of X^^^ restricted to x > coincides with (4). The 
expression of the Levy measure for x < follows easily by symmetry considerations. This proves that 
ijj{z) defined by (1) is the Laplace transform of a (possibly killed) Levy process with the density of the 
Levy measure given by (4). 

The rest of the proof is rather straightforward. Property (v) follows easily from the definition of the 
Laplace exponent (1). The Wiener-Hopf factorization (iv) follows easily by construction: we know that 
both K,{q, z) and k{q^ z) defined by (5) are Laplace transforms of (possibly killed) subordinators, and the 
result follows from identity —ip{z) = ^(g, —z)k{q, z) and the uniqueness of the Wiener-Hopf factorization 
(see [27]). 

Let us prove (ii). The fact that X drifts to +00 when /3 = 1 and /3 > follows from the Wiener-Hopf 
factorization (iv): in this case K{q,0) = 0, therefore the ascending ladder height process drifts to +00, 
while the descending ladder height process is killed at a rate k{q, 0) > 0. The expression for E[Xi] follows 
from (1) using the fact that E[Xi] = ip'{0). Other results in (ii) can be verified in a similar way. 

Let us prove (iii). Formula (1) and the asymptotic expansion for the Gamma function (see formula 
8.328.1 in [24]) imply that 

^^(i^) = 0(1^1^+^) = od^n, z^oo.z e M. (6) 

Applying Proposition 2 from [5] we conclude that X has no Gaussian component, and that when 7+7 < 1 
{l<7 + 7<2} the process has paths of bounded variation and no linear drift {paths of unbounded 
variation}. In the remaining case 7 + 7 = 1 the density of the Levy measure has a singularity of the 
form Cx'"^ + o(x~^) as x — 0+ (see formula 15.3.12 in [1]), which implies that the process has paths of 
unbounded variation. □ 



Remark 1: It is possible to extend the definition of the hypergeometric processes to the case when 
7 G {0, 1} or 7 G {0, 1}. In the case when 7 = 1 and 7 G (0, 1) the process X is spectrally negative. 
This is true due to the formula (4), which shows that 7r(x) = for x > 0; one can also establish this 
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fact using (5), which imhes that z) = 1 — P + z. Similarly, when 7 = 1 and 7 G (0, 1) the process 
X is spectrally positive. When 7 = 7 = 1 the process X is Brownian motion with drift, killed at an 
exponential time. When 7 = {7 = 0} the process —X {X} is a subordinator. 

The three Lamperti-stable processes ^*,^^ and were introduced by Caballero and Chaumont [10] 
by applying the Lamperti transformation (see [30]) to the positive self-similar Markov processes related 
to stable process. In particular, the process ^* is obtained from a stable process started at x > and 
killed on the first exit from the positive half-line, while the process {^^} is obtained from stable process 
conditioned to stay positive {conditioned to hit zero continuously}. We refer to [10, 18, 13] for all the 
details on these processes. 

For our next result we will need the form of the Laplace exponent of the Lamperti-stable process 
We assume that the characteristic exponent "^viz) = — lnE[exp(izYi)] of a stable process Y is defined 
as 

^^(^) =e^z°l{,>o} + e-'?|^|n{,<o}, (7) 

where I7I < a if a G (0, 1) and I7I < 1 if a G (1,2). With this parametrization the density of the Levy 
measure of Y is given by 

7ry(x) = c+a:"^""l{^>o} + c_|x|"^""l{^<o}, 

where 

c+ = r(l + a)^^^^, c_ = r(l + a) ^^"^""^^-^^^ (8) 
TT vr 

and p = FiYi > 0) = (1 — 'j/a)/2. Note that this parameter was corresponds to 1 — p in [10]. According 
to Caballero and Chaumont [10], the Laplace exponent ip*{z) = lnE[exp(2;^J')] of the Lamperti-stable 
process ^* is given by 

^*(^) = ^+~^~ z+ I (e"^ - 1 - 2(e^ - l)l{|e.-i|<i}) e^TTy (e^ - 1) dx - c_a-\ (9) 
1 — Jm.\{Q} 

It is important to note that when a < 1, the Laplace exponent (9) can be rewritten as 

xIj*{z) = [ (e"" - 1) e^TTy (e" - 1) dx - c^a-\ 
Jr\{o} 

so that in this case ^* is a process of bounded variation with no linear drift. 

Theorem 1. Lamperti-stable processes C,'^ can be identified as hypergeometric processes with the 

following sets of parameters 





/3 


7 


(3 


7 


r 


l-a(l-p) 


ap 


l-a(l-p) 


a(l - p) 


e 


1 


ap 


1 


a{l- p) 







ap 





"(1 - P) 
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Proof. For the proof of the result for ^'^ see Proposition 2 in [29]. The result for ^-^ follows from Proposition 
1 in [18]. Thus we only need to prove the result for 

Let us set 7, f3, 7) = (1 — a(l — p), ap, 1 — a{l — p), a{l — p)) and compute the Levy measure of 
the hypergeometric process X. We find that r] = 1 + a, and formula (4) implies that for x > we have 

^ -e-°"-2Fi (1 + ap, 1 + a; 1 + ap; e"") 



r(l + ap)T{—ap) 



r(l + a 



sin(7rap) 



-x\—l—a 



e VTy I e 



TT 



where we have used the reflection formula for the gamma function and the fact that 2-^1(0, 6; a, z) = 
(1 — z)~''. Similarly, we find that 7r{x) = e^nyie^ — 1) for x < 0. We see that the Levy measure of the 
hypergeometric process is the same as the Levy measure of the Lamperti-stable process. We know that 
in the case a < 1 both of these processes have paths of finite variation and no linear drift, this proves 
that X = in the case a < 1. 

In the case a > 1 processes X and ^* have infinite variation, no Gaussian component and identical 
Levy measures. Thus their Laplace exponents may differ only by a linear function. Thus in order to 
establish that the Laplace exponents are equal it is enough to show that E[Xi|(^ > 1] = ]E[^J'|(^ > 1]. 
Using (1) and the reflection formula for the gamma function we find that the Laplace exponent of the 
hypergeometric process X is given by 



= — r(tt — z)T{l + z) sin(7r(z — a(l — p))) 
n 



Therefore we have 



E[Xi|C > 1] = ^'x{0) = r(^) Sm(7ra(l P)) ^^^^^ _ ^^^^^ ^ ^^^^ cos(7ra(l - p)). (10) 



TT 



where "^{z) = r'{z)/r{z) is the digamma function. On the other hand, equation (11) in [10] tells us that 

ln(2) 00 



E[eiic>i] 



a 



+ c+ 



+ c_ 



y 



(gy - 1)1+^ 



-dy 



ln(2) 



[ev - i)i+< 



:iii 



(1 - ev) 



where c+ and c_ are defined in (8). Let us show that the expressions in the right-hand side of (10) and 
(11) are equal. 

First we will deal with the integrals multiplying in (11). We rearrange the terms as follows 



ln{2) 



;i + y-e^)e^ 
(gy - 1)1+° 



■dy 



ye" 



ln(2) 



dy 



[ey 



\1 + Q 



i)ey 



ln(2) 



{ey - 1)1+" 



dy 



(12) 



Performing the change of variables u = exp{y) it is easy to see that the second integral in the right-hand 
side of (12) is equal to 1/(q; — 1). In order to compute the first integral we use integration by parts 



[l + y-ey)ey 
{ey - 1)1+^ 



-dy 



1 

a 



y - e^)d(e^ 



1 

a 



1 - ey 
{ey - ly 



TT 



a sm vra 



(13) 
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where in the last step we have used the following integral formula 

ir(c-f)r(i-c) a 



-du 



b r(i-f) 



< c < 1, 



(14) 



and the reflection formula for the gamma function. The last integral in (11) can also be computed using 
integration by parts 



[1 + y- ey)e^ 
(1 - eyy+'^ 



dy 



[l-y-e-y)e' 



ay 



-dy 



{ev - 1)1+ 



oo 

-- l{l-y- e-^)e("-i)M(e^ - 1)"" 
a J 



oo 

1 /■ (-1 + e-y)e^'^-^^y + {l-y- e-^)(a - l)e^''-^^y 



(15) 



a 



a-2 



a 



{ey - ly 



dy 



{ey - 1) 



a 



ye 



{a-l)y 



a - 2 « - 1 d 

-B(2 — a, 1) — 



a 



a dz 





oo 



{ey - ly 



{ey - 1)" 



dy 



J z=a—l 



= ----(*(!)- ^(2 -«)) 
a a 

where in the last step we have again used (14). Combining (11), (12), (13) and (15) we see that 

1 



1 — a 



a — 1 a sinfTra) 



----{^{l)-^{2~a)) 
a a 



Using the fact that \E'(2 — a) = \E'(a) + vr cot(7ra) + 1/(1 — «), the identity r{z + 1) = zr{z) and addition 
formula for the sine function it is not hard to reduce the above expression to (10), we leave all the details 
to the reader. □ 



3 Mellin transform of the exponential functional 

Let X be the hypergeometric Levy process with parameters {(3, 7, /3, 7) E A and /3 > 0. We assume that 
a > and define the exponential functional as 

C 

I{a,X) = j e-"^'dt, (16) 
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where ( is the hfetime of the process X. Note that the above integral converges with probabihty one: 
according to Proposition 1, either ( is finite (if /3 < 1) or ^ = +00 and the process X drifts to +00 (if 
/3 = 1). 

Our main tool in studying the exponential functional is the Mellin transform, which is defined as 

M{s) = A^(s;a,/3,7,/3,7) = E[I{a,Xr~']. (17) 

From the definition of the Laplace exponent (1) we find that X satisfies Cramer's condition E[exp(— /3Xi)] = 
1, therefore applying Lemma 2 from [39] we conclude that Ai{s) exists for s G (0, 1 + (36). 

In order to describe our main result in this section, we need to define the Barnes double gamma 
function, G{z; r). This function was introduced by Alexeiewsky in 1889 and was extensively studied by 
Barnes in [2] and [3]. The double gamma function is defined as an infinite product in Weierstrass's form 

G(z;T) = -e^^+''^Y\Y\'(l + — ) e"^+2(J+n)^^ |arg(r)|<7r, z e C. (18) 

T J. J. J. J. \ fyij- -\- n J 

m>On>0 ^ ^ 

Here the prime in the second product means that the term corresponding to m = n = is omitted. Note 
that by definition G{z; r) is an entire function in z and has simple zeros on the lattice mr + n, m < 0, 
n < 0. The function G{z; r) can also be expressed as an infinite product of gamma functions, see [2]. An 
integral representation for ln(G(2;; r)) was obtained in [31] and several important asymptotic expansions 
were established in [8]. 

It turns out that it is possible to define constants a = a(r) and b = hij) in a particular way, so that 
we have (7(1; r) = 1 and that G{z\t) satisfies the following three functional identities (see [2] and [25]) 

G{z + 1-t) = r(^)G(2;r), (19) 

G{z + t-t) = (27r)^r^^+^r(2)G(^;r), (20) 
G{z;t) = (2vr)t(W)^-fl+i(i+^)-iG(f;i). (21) 

Next we introduce a function which will play a central role in our study of exponential functionals. 
Everywhere in this paper we will use notation 

^ = i. 

a 

Definition 1. For s & C we define 

M(..) ^ Mis; a, ^ C ^ ' f - ^»;,t;)^; ' . „ . C. ,22) 

G{{1- (3 + 'j)d + s;d) G{/36 + 1 - s; 6) 

where the constant C is such that M{1) = 1. 

Our main result in this section is the following Theorem, which provides an explicit expression for 
the Mellin transform of the exponential functional. 

Theorem 2. Assume that a> 0, {f3,-f,f3,^) G A and /3 > 0. Then M{s) = T{s)M{s) for all s G C. 
Before we are able to prove Theorem 2, we need to establish several auxiliary results. 
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Lemma 1. 

(i) Assume that r > 0. When s — )■ oo m the domain \ arg(s)| < n — e < n we have 



In 



-sln(s) - -(1 + ln(r))s + ^{a - 1 - r) ln(s) + 0(1). (23) 
r r 2r 



G{s; t) 

(iv) When s — oo m the domain < e < arg(s) < tt — e we have 

ln(M(s)) = -(7 + l)s ln(s) + s((l + ln(5))(7 + 7) + 7ri7) + 0(ln(s)). (24) 

(Hi) When s ^ 00 in a vertical strip a < Re(s) < b we have 

\M{s)\ = exp (|(7 - 7)|Im(s)| + 0(ln |Im(s)|)) . (25) 

Proof. Part (i) follows from the asymptotic expansion for G{z;t), given in formula (4.5) in [8], while 
parts (ii) and (iii) are simple corollaries of (i) and (22). □ 

The following notation will be used extensively in this paper: if X is a hypergeometric process with 
parameters (/3,7,/3,7), then X denotes the hypergeometric process with parameters {6(3, S'j, 6(3,6^), 
provided that this parameter set is admissible. In particular, the Laplace exponent of X is given by 

r{i-6p-z)r{6p + z) 

Lemma 2. 

(i) M{s) is a real meromorphic function which has zeros 

- {1 - /3)S - mS - n, 1 + + j)S + mS + n\ , (27) 

y m,n>0 



and poles 



-{1 - (3 + j)6 -m6 -n, z:^^^ = 1 + (36 + m6 + n} . (28) 

m,n>0 



All zeros/poles are simple if a ^ Q. 
(ii) M{s) satisfies the following functional identities 



M{s + 1) = -^—M{s), (29) 

-<5(7+7) 

Mis + 5) = -7—- -M(.), (30) 

ip{l — d — s) 

M{s;a,(3,-fJ,^) = a^^-'^^^+''^M{l - a + as;6,6(3,6-f,6P,6^). (31) 
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Proof. The proof of (i) follows from the definition of the double gamma function (18), while functional 
identity (29) is a simple corollary of the functional identity for the double gamma function (19). 
Let us prove (31). We use (21) and find that for all s and x 

G{s + x; 5) ^ ^asx^ G{a{s + xy,a) 
G (s; 6) G{as; a) 

where C = C{x) depends only on x. The above identity implies that 

G{{1- P)6 + s;6) Gi0 + ^)6 + l- s;6) 
GUI- P + -f)6 + s-6) G0S+l-s;6) 

_ ^(1 - P + as; a) G0 + 7 + a - as; a) 

(7(1 - /3 + 7 + as; a) G0 + a-as]a) 

where C does not depend on s. It is easy to see that (31) follows from the above identity and (22). 
The functional identity (30) follows from (31) and (29), we leave all the details to the reader. □ 



The following proposition will be central in the proof of Theorem 2. It allows us to identify explicitly 
the Mellin transform of the exponential functional. This result is also applicable to some other Levy 
process, including Brownian motion with drift and more generally, processes with hyperexponential or 
phase-type jumps, therefore it is of independent interest. 

First let us present the main ingredients. Let y be a (possibly killed) Levy process started from 
zero, and let iPy{z) = lnE[exp(zFi)] denote its Laplace exponent. In the case when V'y(O) = (the 
process is not killed) we will also assume that E[Yi] > 0, so that Y drifts to +00. As usual we define 
the exponential functional / = exp{—Yt)dt (where ( is the lifetime of Y) and the Mellin transform 
A^y(s) =E[/^-i]. 

Proposition 2. (Verification result) Assume that Cramer's condition is satisfied: there exists zo < 
such that ipvi^) is finite for all z G {zq,0) and ipyi—d) = for some 9 G (0, — -^o)- U f{s) satisfies the 
following three properties 

(i) /(s) is analytic and zero-free in the strip Re(s) G (0, 1 + 6'), 

(ii) /(I) = 1 and f{s + 1) = -sf{s)/iljY{-s) for all s G (0, 6), 

(ill) |/(s)|-i = o(exp(27r|Im(s)|)) as Im(s) 00, Re(s) G (0, 1 + 9), 

then My{s) = /(s) for Re(s) e{0,l + 9). 

Proof. The Cramer's condition and Lemma 2 in [39] imply that A^y(s) can be extended to an analytic 
function in the strip Re(s) G (0, 1 + 9). In the case when iPy{0) = {^/^y(0) < 0} we use Lemma 2.1 in 
[32] {Proposition 3.1 from [15]} to conclude that A^y(s) satisfies the functional identity 

My{s + 1) = -^-My{s) (32) 

for all s G (0,6*). Since /(s) satisfies the same functional identity we conclude that the function F{s) = 
■M.y{s) / f{s) satisfies F(s + 1) = F{s) for all s G (0,6*). Using the assumption that /(s) is analytic and 
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zero-free we conclude that F{s) is an analytic function in the strip Re(s) G (0, 1 + 6*). Since F{s) is also 
periodic with period equal to one, it can be extended to an analytic and periodic function in the entire 
complex plane. 

Our goal now is to prove that function F{s) is constant. Since F{s) is analytic and periodic in the 
entire complex plane, it can be represented as a Fourier series 

F{s) = Co + Fi(exp(27ris)) + F2(exp(-27rzs)), (33) 

where Fi[z) are entire functions 

Fi(^) = 5^c„z", F2(z) = ^c_„z" 

n>l n>l 

Due to the inequality |A^y(s)| < A^y(Re(s)) and assumption (iii) we conclude that F{s) = o(exp(27r|Im(s)|)) 
as Im(s) — )■ oo. As Im(s) — )■ — oo this estimate and equation (33) imply that Fi(z) = o{z) as 2; — )■ oo, 
and using Cauchy's estimates (Proposition 2.14 in [20]) we conclude that -^1(2;) = 0. Similarly we find 
that -^2(2) = 0, thus F{s) is a constant, and the value of this constant is one, since A^y(l) = /(I) = 1.0 



We would like to stress that Proposition 2 is an important result of indepent interest. We know that 
if Cramer's condition is satisfied then the Mellin transform Aiyis) satisfies the functional identity (32), 
however it is clear that there are infinitely many functions which satisfy the same functional identity. 
Proposition 2 tells us that if we have found such a function /(s), which satisfies (32), and if we can verify 
the two conditions about the zeros of this function and its asymptotic behaviour, then we can in fact 
uniquely identify A^y(s) = f{s). In particular, this proposition can be used to provide a very simple 
and short proof of the well-known result on exponential functional of Brownian motion with drift and 
of the recent results on exponential functionals of processes with double-sided hyper-exponential jumps 
(see [14]). 

Proof of Theorem 2: First of all, we check that the Cramer's condition is satisfied with Q = pS. Let 
f{s) = r(s)M(s), where M{s) is defined by (22). From Lemma 2(i) we know that f{s) is analytic and 
zero-free in the strip Re(s) G (0, 1 -|- (36). By construction we have /(I) = 1, and from the formula (29) 
we find that /(s) satisfies f{s + 1) = —sf{s)/ip{—as) for s G (0,(36). Next, Lemma l(iii) and Stirling's 
asymptotic formula for the gamma function (see formula 8.327.3 in [24]) 

ln(r(s)) = s ln(s) - s + 0(ln(s)) (34) 

imply that as s — 00 in the vertical strip Re(s) G (0, 1 + P6) we have 

\f{s)r' = exp (I (1 - 7 + 7) |Im(s)| + o(Im(s))) = a (exp(7r|Im(s)|)) , 

where in the last step we have also used the fact that both 7 and 7 belong to the interval (0, 1). 

We see that function f{s) satisfies all conditions of Proposition 2, thus we can conclude that Ai{s) = 
f{s). □ 



Corollary 1. Assume that a > 0, (3 > and that both sets of parameters {(3, 7, (3, 7) and {6(3, ^7, 6(3, 6^) 
belong to the admissible set A. Then we have the following identity in distribution 

X I{a; X) = x ei x I{6; X)" (35) 
where ei ~ Exp{l) and all random variables are assumed to be independent. 
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Proof. Rewrite (31) as 

r(l - a + as)M{s; a, /3, 7, /3, 7) = a^'-^^^"'+''^T{s)M{l - a + as; S, S(3, S-f, 5/3, ^7) 

and use the following facts: (i) r(s) = E[e^~^]; (ii) if f{s) is the Mellin transform of a random variable 
1^ then /(I — a + as) is the Mellin transform of the random variable (iii) Mellin transform of the 
product of independent random variables is the product of their Mellin transforms. □ 



4 Density of the exponential functional 

In this section we will study the density of the exponential functional, defined as 

d 



p{x) 



dx 



P(/(a,X) < a;), a; > 0. 



As in the previous section, X is a hypergeometric process with parameters 7, /3, 7) G ^ and f3 > 0. The 
main results of this section are the convergent series representations and complete asymptotic expansions 
of this function as x — )■ 0^ or x — )■ +00. 

Let us define the following three sets of parameters which will be used extensively later. In the 
following definition (and everywhere else in this paper) ipi') ^ind M(-) denote the functions which 
were defined in (1), (26) and (22); the sequences 2;"^ and z:^^ represent the poles of M(-) and were 
defined in (28); the constant rj is defined by (3). 



Definition 2. Define the coefficients {a„,}„>o as 

1 " 

On = ^JJ^(«i), n > 0. 



(36) 



j=0 



The coefficients {bm,n}m,n>o are defined recursively 



^0,0 



6 



/3 + 7)5) 



r(r?-7)r(-7) 



M(l- (1-/3 + 7)5), 



'm,n-l, m > 0, n > 1, 
T(z~ ) 



(37) 



T(z. 



-'m—l,ni 



m > 1, n > 0. 



Similarly, {cm,n}m,n>o are defined recursively 

T(l + /35)r(l-/3 + /3) 



Co,0 



r(r/-7)r(7) 



m—l,n) 



M05) 



7 

rc™,„-i, ,m>0, n>l. 



^(-«Vn-l) 



(38) 



a 



-5(7+7) 



Tizl 



V^(l - 5 - J T{zt.,,n) 



Cm_i n, m > 1, n > 0. 
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Note that if /3 = 1 we have ip{0) = 0, which imphes that a„ = for all n > 0. 
Proposition 3. Assume that a ^ Q. For all m,n > we have 

Res(A^(s) : s = —n) = a„, if /3 < 1, 
Res{M{s) : s = z- J = K 
Res(A^(s) : s = z^n) = —c„ 



''m,ni 



Proof. Let us prove that the residue of A^(s) at s = z^^ is equal to bm,n- First, we use Theorem 2 and 
rearrange the terms in the functional identity (29) to find that 

, ^ , , 6M{s + 1) r(2 - /3 + 7 + as)T0 + 7 - as) 

M.(s) — 



s + (1-/3 + 7)5 T{1- (3 + as)T0-as) 
The above identity and definition (37) imply that as s — )■ — (1 — /3 + 7)5 

which means that the residue of A^(s) at z^q = —(1 — /3 + 7)^ is equal to 6o,o- 

Next, let us prove that the residues satisfy the second recursive identity in (37). We rewrite (29) as 

= _^iz^7U(s + l). (39) 

s 

We know that (s) has a simple pole at s = 2;" ^ while (s + 1) has a simple pole at „ + 1 = z^_i 
One can also check that function ipl—as) is analytic at s = z"^ for m > 1. Therefore we have as 

M{s) = Res{M{s) : s = z-J _^ „ +0(1), 
Mis + 1) = Res{Mis):s = z-_,J-^ + Oil), 

: — = — + 0{s-z^j 



which, together with (39) imply that 

?/'(— az~^) 

Res(A^(s) : s = 2„„) = — '— x Res(Al(s) : s = z^_^^J. 

The proof of all remaining cases is very similar and we leave the details to the reader. □ 



Proposition 3 immediately gives us a complete asymptotic expansion oip{x) as x — )■ O"*" and x — )■ +00, 
presented in the next Theorem. 
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Theorem 3. Assume that a ^ Q. Then 

n>0 

' EE' 



p[x) 



m>0 n>0 

{m+f^)S-n-l 



X — +00. 



(40) 
(41) 



m>0 n>0 

Proof. The starting point of the proof is the expression of p{x) as the inverse Melhn transform 



p{x) = — [ M{s)x Ms, X > 0. 
2m J 



(42) 



l+iB 



Due to (25), Theorem 2 and Stirhng's formula (34) we know that |A^(x + iu)\ decreases exponentially 
as u — 00 (uniformly in x in any finite interval), therefore the integral in the right-hand side of (42) 
converges absolutely and p{x) is a smooth function for x > 0. Assume that c < satisfies c ^ z:^^ and 
c 7^ —n for all m, n. Shifting the contour of integration 1 + iM 1— )■ c + iM and taking into account the 
residues at the poles s = 2;" „, we find that 



p{x) = y^Res(A^(s) : s = z" ,,) x x"''"'-" + ^ Res(A^(s) : s 

0<n<\c\ 
1 



-n) X X 



(43) 



+ 



27ri 



M(s)x Ms. 



c+iR 



where the first summation is over all m > 0,n > 0, such that 2;„„ > c. Next, we perform a change of 
variables s = c + iw and obtain the following estimate 



M(s)x Ms 



X 



A^(s)x~™dM 



< x' 



J |A^(c + iM)|dM = 0{x' 



which proves (40). The proof of (41) is identical, except that we have to shift the contour of integration 
in the opposite direction. □ 



It turns out that for almost all parameters a the asymptotic series (40) and (41) converge to p{x) for 
all X > 0. In order to state this result, we need to define the following set of real numbers. 

Definition 3. Let C be the set of real irrational numbers x, for which there exists a constant b > 1 such 
that the inequality 



p 

X 



1 



(44) 



is satisfied for infinitely many integers p and q. 



This set was introduced in [25] in the connection with the distribution of the supremum of a stable 
process and it was later studied in [23]. It was proved in [23] that x G £ if and only if 



lim 

g— >+oo 



In(gx) 



0. 



(45) 
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where (x) = min{|x — i\ : i G Z}. There also exists a characterization of elements of C in terms of 
their continued fraction expansion (see Proposition 1 in [23]). It is known that £ is a proper subset of 
Liouville numbers, that it is dense in M and that it has Lebesgue measure zero. The Hausdorff dimension 
of C is also zero. This set is closed under addition/multiplication by rational numbers. It is also known 
that X G £ if and only if G C. See [23] for proofs of these results and for some further references. 
The following Theorem is our main result in this section. 

Theorem 4. Assume that a ^ £ U Q. Then for all x > 

E anx^ + E E if 7 + 7 < 1, 

^(^\ ^ ; «>0 m>On>q /^gN 

' ^ E E S if 7 + 7 > 1. ^ ^ 

m>0 n>0 

First let us establish the following technical result, which gives us a formula for M{s) similar to the 
reflection formula for the Gamma function. 

Lemma 3. Define 

Ck = l-{l-l3 + ^)6-6/2-k. (47) 

Then for all u E C 

k—l 

. ./ N , ^L.,., , N -i-T cos(7ra(? — iw + 7^)) F(—iu) , 

M{ck + iu) = {-l)'M{ck + k + m)xl[ . X \ / 48 

~"_o cos(7ra(j — m)) F[—m + k) 

where we have defined 

, . . , . GiU + w;5) G((r7-7 + i)5 + w;5) , , 

Proof. Let us set s = Ck + iu. Then iterating the functional identity (29) k times we obtain 

^ r(s + A:)J-jL r(l-/5 + «s + «j) T{f3-as-aj) ^ ' 



We use identity 



T{s + k) ' ' r(l-s) 
and the reflection formula for the gamma function and rewrite (50) as 



, n(l-fl + l)l> + SI2-m) 
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Next we change the index j ^ k — 1 — j and use reflection formula for the gamma function to obtain 



^ V{-a{k - 1 - j) - I + iau) _ t-^ cos(7ra(j -'m + 7^)) r(aj + 7 + | - iau) 
V{-a{k - l-j)-7-| + iau) cos(7rQ;(j - iu)) r(aj + f - mu) 

The proof of (48) follows from (51), (53) and the following identity 

n''"^ . , X G{5z + k-5) 

which is obtained by iterating formula (19) k times. □ 



Proof of Theorem 4: We will use a similar technique as in the proof of Theorem 2 in [23]. Let us 
define 5 = 1—7 — 7 assume that B > 0. We start with (43) and set c = c^, where is defined 
by (47). Note that A4{s) does not have singularities on the vertical line Ck +iM; if this was not the case 
then Cfc would coincide with one of the poles ^~ „, which would imply that a is rational. 

We define 



h{x,k) = x-'^Re 



M{ck + m)x '""du 



/2(x, k) = X "''Re 



M{ck + m)x '"dw 



(54) 



It is clear that the integral in the right-hand side of (43) is equal to 2{Ii{x, k) + l2{x, k)). Our goal is to 
prove that Ij{x, A;) — )■ as /c — )■ +00 for all x > 0. 

First, let us deal with l2{x, k). Using Theorem 2, formula (24) and Stirling's asymptotic formula (34) 
we find that there exists a constant Ci > such that for all s in the domain 

V = {s eC : \s\ > 2 and - < arg(s) < — } 

8 8 

we have an upper bound 

\M{s)\ < \sf' exp (Re [Bs\n{s) + s((l + ln(5))(l - B) + - 1)]) . 

From now on we will denote u = Im(s). Computing the real part in the above expression we obtain that 
for all s G P 

\Mis)\ < |s|^iexp(5Re(s)ln|s|-(7r7 + 5arg(s))M + Re(s)((l + ln(5))(l-5)-l)) (55) 

Next we check that for all k large enough, conditions u > k and s = Ck + lu imply s G so the 
integrand in formula (54) defining /2(x, k) can be bounded from above as given in (55). Let us restrict 
s to the line of integration L2 = {s = Ck + lu : u > k} and simplify this upper bound. Basically, 
we want to isolate a term of the form exp(— 5/;; ln(/c)) and show that everything else does not grow 
faster than an exponential function of k. First of all, when s G L2 we have Re(s) = < 1 — /c and 
|s| > |Re(s)| = \ck\ > A; — 1, therefore 



exp (5Re(s) In |s|) < exp(-5(A; - 1) ln(A; - 1)). 
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Next, for s e L2 and k sufficiently large it is true that —2k < Re(s) < 1 — k, which implies that there 
exists a constant C2 > such that 

exp (Re(s)((l + ln(5))(l - B) - 1)) < Cl 

Finally, for s E L2 we have arg(s) > 0, which together with the assumption B > shows that for all 
u>0 

exp (— (vr7 + B a,Tg{s))u) < exp{—7i^u) . 
Combining the above three estimates with (55) and (54) we see that 



\l2{x,k)\ < x-'" J \M{ck + m)\du 

k 



(56) 



The right-hand side of the above inequality converges to as /c — )■ +00, therefore l2{x,k) — )■ as 
k — )■ +00. 

Now we will deal with /i(x, k). Our ffist goal is to find an upper bound for the product of trigonometric 
functions in (48). We will follow the proof of Theorem 2 in [23]: we use the trigonometric identities 

cos(x + iy) = cos(x) cosh(?/) — i sin(x) sinh(y), 
I cos(x + iy)|^ = cosh(2/)^ — sin(a;)^, 

which imply that | cos(x)| cosh(i/) < | cos(x + i?/)| < cosh(?/), therefore 



cos(a + iy) 



cos(6 + iy) 



< 



cos 



Applying the above estimate and Lemma 1 from [23] we conclude that for a ^ £ U Q and for k large 
enough 



fc-i 

n 

j=0 



cos(7ra(j — iu + 7^)) 



cos(7ra(j — iu)) 



k-l 



< JJ|sec(7r«j)l =2'=+°('=) < 3^ 

j=0 



Using (54), (50) and the above inequality we conclude that for all k large enough 

k 



/i(x, k)\ < x-'"'3'' I \M{ck + k + in) I x 



F{-iu) 



F{-iu + k) 



du. 



(57) 



Now our goal is to prove that F{iu) / F{—iu + k) converges to zero faster than any exponential function 
of as A; — 7- +00. We use (49), Stirling's formula (34) and the asymptotic expansion (23) to conclude 
that when w — )■ 00 in the domain | arg(w)| < 37r/4 we have 



ln{F{w)) = Bwln{w) + 0{w) 
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This asymptotic result implies that there exists a constant C3 > such that for k large enough and for 
all V e [0, 1] 



F{-ikv) 



F{-ikv + k) 



< C3 exp (5Re [{-ikv) \\i{-ikv) - k{l - iv) ln(A;(l - iv 

= C3 exp (^B -k \n{kVl + f^) + kv ^arctan 

< C^exp{-Bk\n{k)) 



[v] 



where in the last step we have used the fact that arctan(f) < tx/2 and -B > 0. Thus we see that for all 
k large enough we have 



max 

0<u<k 



F{-m) 



F{-m + k) 



max 

0<v<l 



F{-ikv) 



F{-ikv + k) 



< C^exp{-Bk\n{k)) 



Combining the above estimate with (57) we obtain that for all k large enough 



Ji(x, k)\ < a;^"*3^Cfe^^^'"('=) / \M{1 - (1 - /3 + 7)5 - 6/2 + m)\du 



and we see that fc) — )■ as — +00. Thus when 7 + 7 < 1 the first series in (46) converges to p{x) 
for all X > 0. 

When 7 + 7 > 1 the proof is very similar, except that one has to shift the contour of integration in 
the opposite direction. □ 



5 Applications 

In this section we will present several applications of the above results on exponential functionals of the 
hypergeometric Levy processes. In particular, we will study various functionals related to strictly stable 
Levy processes. Our main tool will be the Lamperti transformation, which links a positive self-similar 
Markov process (pssMp) to an associated Levy process. By studying this associated Levy process and 
its exponential functional we can obtain many interesting results about the original self-similar Markov 
process. 

We assume that the reader is familiar with the Lamperti transformation [30] and Lamperti-stable 
processes (see [10, 18, 13, 12]). 



5.1 Extrema of stable processes 

As our first example, we will present a new proof of some the results on the distribution of extrema of 
general stable processes, which were obtained in [25] and [23]. Let us define 5*1 = sup{Yi : < t < 1}, 
where F is a stable process with parameters (a,p), whose characteristic exponent is given by (7). The 
infinite series representation for the density of Si (valid for almost all values of a) was derived recently 
in [23]. This formula was obtained using an explicit expression for the Mellin transform of Si, which was 
found in [25] via the Wiener- Hopf factorization of stable processes. Our goal in this section is to present 
a more direct route which leads to the density of Si. 
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Proposition 4. We have the following identity in distribution 



^r"^J(a,0, (58) 
where ^ is a hypergeometric Levy process with parameters [ap, ap, ap, a{l — p)). 

Proof. Some of the arguments in our proof have been used before in [18], nevertheless, we provide them 
here for completeness. Let us consider the dual process Y = —Y. Note that the process Y is itself a 
stable process with parameters (a, 1 — p). Let Z be the process killed at the first exit from the positive 
half-line 

Yt, ift<f, 
A, if t>f, 

where T = inf{t : Yt < 0} and A is the cemetery state. Let ^ be the process associated with Z by the 
Lamperti transformation. We use Theorem 1 and find that ^ is a hypergeometric process with parameters 
(1 — ap, a{l — p), 1 — ap, ap). 

Next, from the Lamperti transformation we find that T = /(a,— ^). This statement is equivalent 

to T = I{a,$,), where = — is a hypergeometric process with parameters {ap,ap,ap,a{l — p)) (see 
Proposition l(v)). The rest of the proof follows from the chain of identities 

Pi(f > t) = Pi( inf Ys>0)= Po( inf > -1) 

0<s<t 0<s<t 

= Po( sup < 1) = Po(t'/" sup y; < 1) = Po(5r° > t). 

0<s<t 0<s<l 

□ 



Proposition 4 implies that the density of 5*1 can be represented in terms of the density of I{a,C,) 
(denoted by p{x)) as follows 

-^nSi <x) = ax-^-^pix-"), x>0. 
ax 

Using this expression and Theorem 3 {Theorem 4} we recover the asymptotic expansions that appears 
in Theorem 9 in [25] {the convergent series representations given in Theorem 2 in [23]}. 

Remark 2: Note that the convergent and asymptotic series given in Theorems 4 and 3 have identical 
form. Combining these two results we obtain the following picture: depending on whether 7 + 7 < 1 
or 7 + 7 > 1 one of the series in (46) converges for all x G (0,oo), and this convergent series is also 
asymptotic at one of the boundaries of this interval (at or +00). At the same time the second series in 
(46) provides an asymptotic series at the other boundary. Therefore in what follows we will only present 
the statements about the convergent series representation, with the understanding that the asymptotic 
expansions are implicitly imbedded into these results. 
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5.2 Entrance laws 



In this section we will obtain some explicit results on the entrance law of a stable process conditioned to 
stay positive and the entrance law of the excursion measure of the stable process reflected at its infimum. 
We will use the following result from [6]: if a non-arithmetic Levy process X satisfies E[|Xi|] < oo and 
E[Xi] > 0, then as x — >■ 0+ its corresponding pssMp (Z, Q^.) in the Lamperti representation converges 
weakly (in the sense of finite-dimensional distributions) towards a nondegenerate probability law (Z, Qo). 
Under these conditions, the entrance law under Qo is described as follows: for every t > and every 
measurable function / : M+ — M+, 



(59) 



where / = /(a,X) = Jq°° exp{— aXs} ds. Necessary and sufficient conditions for the weak convergence 
of (Z, Qj,) on the Skorokhod space were given in [11]. 

As before, we consider a stable process {Y^F^) and denote by {Y,¥\) the stable process conditioned 
to stay positive (see [16, 17] for a proper definition). It is known that (F, P|.) is a pssMp with index 
a. According to Corollary 2 in [10] (and Theorem 1 above) its associated Levy process is given by 
which is a hypergeometric Levy process with parameters (1, ap, 1, a(l — p)). From Proposition l(ii) we 
see that IE[^[] > 0, and the fact that hypergeometric processes have Levy measures with exponentially 
decaying tails implies that IE[|^J|] < oo, thus we have the weak convergence of (y, P|.) as x ^ O"*". We 
denote the limiting law by P^. Note also that in this particular case, the weak convergence of (F, P]J) has 
been proved in a direct way in [16]. 

Everywhere in this section the coefficients {&m,n}m,n>o and {cm,n}m,n>o are defined as in Definition 2 
with parameters /3 = /3 = 1, 7 = ap and 7 = a(l — p). 

Proposition 5. Let p\ he the density of the entrance law of (F, P'''). If a ^ CUQ then for x > we 
have 



-i-p 



E E ^m,nX-"^/"-", if a < 1 



Pli^) — { ^i/a m/a+n -f ^ i ^^^^ 

ar{ap)r{l+ail-p)) Z^ Z^^^^nX , It « > i. 

m>0 n>0 

Proof. Let p{x) be the density of /(a,^^). From (59) we find that the density of the entrance law of 
(Y,F'^) is given by 

pi(^) = 7-x~^p(x~^). 

«E[el] 

Using Proposition l(ii) we obtain K[^l] = T{ap)T{l + a{l — p)). In order to finish the proof we only 
need to apply results of Theorem 46. □ 



Proposition 6. Let qt he the density of the entrance law of the excursion measure of the reflected process 
{Y — H, P), where Y_^ = info<s<t Yg. If a ^ CU Q then for x > we have 



-l-p-Q(l-p) 



V V 6 ^~m/a-n if q, < 1 

Qlix) = ^ ^i/c»-a(i-p) ^ m/a+n ■£ ^ i ^^^^ 

EE c^^nx""'^ , if a > 1. 



ar{ap)T{l+a(l~p)) 

m>0 n>0 
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Proof. The proof follows from Proposition 5 and Theorem 3 in [16], which gives us the identity qi = 



5.3 The distribution of the hfetime of a stable process conditioned to hit 
zero continuously 

Our third application deals with stable processes conditioned to hit zero continuously. Processes in this 
class are defined as Doob's /i-transform with respect to the function h{x) = x°'^^~''^~^, which is excessive 
for the killed process {Ytl{t<T},^x)- Its law F^., which is defined on each a-field J-f by 



dP„. 



ya{l-p)-l 

Mt<T}, (62) 



2»q(1-P)-1 



is that of a pssMp that hits zero in a continuous way. According to Corollary 3 in [10] (and Theorem 1 
above), its associated Levy process in the Lamperti representation is the hypergeometric Levy process 
with parameters (0,ap, 0,a(l — p)). In particular from Proposition l(ii) we find that the process 
drifts to — oo. Let T-^ be the life-time of the stable process conditioned to hit zero continuously. 
From the Lamperti transformation it follows that = I{a, —^^), therefore = /(a,^), where ^ is a 
hypergeometric process with parameters (l,a(l — p),l,ap). Let us define the coefficients {bm,n}m,,n>o 
and {cm,n}m,n>o Definition 2 with /3 = /3 = 1, 7 = a{l—p) and 7 = ap. Then according to Theorem 
46, if o; ^ £ U Q and a; > we have 

. f t'-" E E fem,nt™/"+", if « < 1, 

—W^(T^ <t)=< rn>On>0 

dt I r^-i/" E E W"'"/""", if«>i. 

^ m>On>0 

5.4 Distribution of some functionals related to the radial part of a sym- 
metric stable process 

Our last application deals with the radial part of symmetric stable processes in M*^. Let Y = (Yt,t > 0) 
be a symmetric stable Levy process of index a G (0, 2) in {d > 1), defined by 

for alH > and A G M.'^. Here Fy denotes the law of the process Y started from 

y e R'^, II ■ II is the norm 

in M.'^ and (■, ■) is the Euclidean inner product. 

According to Caballero et al. [12], when a < d the radial process Rt = |||Yt|| is a transient positive 
self-similar Markov process with index a and infinite lifetime. From Theorem 7 in [12] we find that the 
Laplace exponent of its associated Levy process C, is given by 

r((a-z)/2) T{iz + d)/2) 
T{-z/2) T{{z + d - a)/2) 

This shows that the process X := 2^ is a hypergeometric Levy process with parameters (1,0/2, {d — 
a)/2,a/2). From Proposition l(ii) we see that ^ drifts to +00 and that 

EI6I = f [Iff . (64) 
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Let Pj: denotes the law of the process R starting from x > 0. The process ^ satisfies IE[^i] > and 
E[|^i| < oo, and according to [6, 11] we have the weak convergence of {R, Px) towards {R, Pq) as x — j- O"*". 

Proposition 7. Let pt be the density of the entrance law of {R,Pq). If a < d then for x > we have 

;^Esin(^)r(l + f)r(i^)^px-, ifa<l, 
P'^''^=< 4,zi±!^% £(f^Mr w« .r^. (65) 



. 2 ; m>0 



Proof. We follow the same approach as in the proof of Proposition 5. Using (59) we conclude that 
pi{x) = (2/q;E[,^i])x~^p(x"^), where p{x) is the density of = I{a/2,X). In the case a ^ C U Q 

the infinite series representation for p{x) and the expression for the coefficients follow from Theorem 46 
and Definition 2 with /3 = 1, 7 = a/2, P = {d — a)/2 and 7 = a/2. Note that both infinite series in (65) 
converge for all a, thus we can remove the assumption a ^ £ U Q. 

There is also a simpler way to derive an explicit series representation for p{x). Theorem 2 and 
Definition 1 tell us that the Mellin transform of I{a/2,X) is given by 

E|/(a/2. AT-'l = r(s)M " 1, " 1^ Cr(.)-;M_^II^^^. (66) 



2- •2- 2 '2J ^ 'G(l + .;|) G (| - »; J) 

where C is a constant such that the right-hand side of the above identity is equal to one for s = 1. Using 
the functional identity (19) for the double gamma function we simplify the right-hand side in (66) and 
obtain 

E[/(a/2,X)^-] = f^r(.)^i^. (67) 

We see that E[/(q!/2, X)'^~^] has simple poles at points {— n}„>i and {{d + 2m) / a}rn>o- The residues at 
these points (which give us the coefficients in (65)) can be easily found using the fact that the residue of 
r(s) at s = —n is equal to (— 1)"'/?t,!. □ 



Next, we will study the last passage time of (Y,Po) from the sphere in of radius r, i.e. 

Ur = sup{s > : ||YJ < r}. 

Due to the self-similarity property of Y we find that random variables Ur satisfy the scaling property 
b"'Ur = Ubr (valid for any b,r > 0), thus it is sufficient to consider the case of of f/2- 



Proposition 8. If a < d then 



d 



^Esin(^)r(^)r(^)(^r, ifa<i, 



Proof. Let us define the last passage time of R as Lr = sup{s > : ||-Rs|| < r}. It is clear that U2 = Li. 
According to Proposition 1 in [19], the random variable Li has the same law as Q°'I[a, C,), where is the 
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Levy process associated to R by the Lamperti transformation and Q is independent of ^. Lemma 1 in 

[19] tells us that Q = e~^^, where U and Z are independent random variables, such that U is uniformly 
distributed over [0, 1] and the law of Z is given by 



oo 



E[H: 

u 

where H = {Ht, t > 0) denotes the ascending ladder height of ^ and u its Levy measure. From the proof 
of Lemma 3 in [12] we know that H has no linear drift and that its Levy measure is given by (up to a 
multiplicative constant) 

Using the above expression, integration by parts, the integral identity (14) and reflection formula for the 
gamma function we obtain 



oo 

E[Hi] = j xi/(dx) 



TT 



a sinful 



Let us find the Mellin transform of Li. We use the independence of U and Z and obtain 



E[^"] = E [e-'"^l = -E 

s 



1 - e~'^A 1 /■ (1 - e-'^)e^^ , 

z - 



sE[i/i] J (e2- - l)i+"/2 



^ /(I - e-)d [(e^^ - 1)""/^] 



asE^ 



oo 



1 [ e-^^ _ 1 r(f(i + s))r(i-f) 



oElHA J (e2- - l)"/2 aE\Hi] F (l + ^) 

\ I J 

where we have again applied integration by parts and integral identity (14). Combining the above two 
expressions with (67) we have 



Ff 1 F ( ^ 



F(^)r(i-^ 

The function in the right-hand side of the above equation has simple poles at points {—n}n>Q and 
{{d + 2m) / a}rn>o- We express the density of Li as the inverse Mellin transform, compute the residues 
and use similar technique as in the proof of Theorem 46 to obtain the series representation (68). □ 
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